RECURSIVE DISCRETE FOURIER TRANSFORMATION METHOD AND 
RECURSIVE INVERSE DISCRETE FOURIER TRANSFORMATION METHOD 



5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an arithmetic operation 
method for Fourier transformation or inverse Fourier 
10 transformation and more particularly to an arithmetic 
operation method for carrying out those transformations with 
a simple arithmetic operation processing for obtaining a result 
of the arithmetic operation by those transformations in a short 
time . 

15 

2 . Description of the Related Art 

Conventionally, the Fourier transformation method used 
in frequency analysis and the like of time- series data stream 
has been employed in high compression encoding of audio signal 

20 and moving picture signal as well as spectrum analysis of 
signals in audio signal processing field, moving picture data 
processing field for medical system and the like. Further, 
inverse Fourier transformation method and Fourier 
transformation method, which operate in complementary way to 

25 each other, have been widely used as modulation technology 
and demodulation technology in communication field. 

According to the Fourier transformation method, data 
stream sampled as digital amount is treated as N groups (N 
is an integer, for example, 1024) and a time interval, in which 

30 the N-data stream exist, is assumed as a window period while 
that window period is assumed to be basic frequency. Then, 
signal components of data stream existing in the window period 
are obtained as real number component and imaginary number 
component of a harmonic signal of the basic frequency. 

35 The inverse Fourier transformation method is a 

transformation method, which is complementary to the Fourier 
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transformation method. 

Then, data stream supplied to be subjected to the Fourier 
transformation is discrete data stream sampled at every 
predetermined period and a method for conducting the Fourier 
5 transformation on the discrete data is called discrete Fourier 
transformation (DFT) . The discrete Fourier transformation 
technology is employed as an analysis technology for control 
technology, in which, for example, a state of production 
process is obtained with discrete data, the quality of the 

10 process is kept optimum by analyzing the obtained data and 
the rate of acceptable products in manufactured products is 
improved. The inverse Fourier transformation method for 
discrete data is called inverse discrete Fourier 
transformation (IDFT), which is employed for modulation of 

15 digital ground broadcasting , scheduled to be started from 2003 
as well as generation of sound source in audio field, while 
its application field has been widening every year. 

The discrete Fourier transformation technology and the 
inverse discrete Fourier transformation technology, which are 

20 employed in the above manner, are signal conversion 
technologies, which operate in complementary manner. In the 
discrete Fourier transformation of the former, supplied 
signals are sampled at every predetermined time interval and 
a voltage value obtained by the sampling is obtained as sampled 

25 data. When the data stream, which is a group of the obtained 
data , is N data obtained in a predetermined time t,x(t),x(t+l), 

x (t + 2) , x(t+N-2) , x(t+N-l) , avaluex(k, t) of the discrete 

Fourier transformation obtained with respect to the N data 
is defined in a following equation. 



30 
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for k = 0,1, ...,N-1 

As evident from this equation, the Fourier 
5 transformation is carried out on supplied data stream just 
such that its inherent basic function is convoluted in every 
point to be obtained and an operation for convoluting the basic 
function is conducted by a plurality of multiplication 
processings . 

10 When the multiplication processing is carried out with 

a special multiplying circuit or a digital signal processor 
(DSP), it has been known that a load upon hardware used for 
arithmetic operation in the multiplying circuit or the DSP 
is increased tremendously. 

15 As for the load on this multiplication processing, a 

transformation equation expressed with the aforementioned 
equation (1) needs multiplications of 4 X N 2 times, so that 
in case where N is 1024, for example, the multiplication 
frequency is about 4,200,000 times. Therefore, the load of 

20 the arithmetic operation is very large and further, if the 
number N of points to be handled as the data stream is increased, 
the multiplication processing frequency is increased by second 
power, which is unfavorable in terms of actual field. 

Thus, as the discrete Fourier transformation equation, 

25 fast Fourier transformation (FFT) , whose operation efficiency 
is raised by deforming the matrix by paying attention to the 
regularity that the basic function is composed of a cyclic 
function, is employed actually. 

An arithmetic operation method called butterfly 

30 operation is employed for the FFT. The butterfly operation 
is so constructed that a simple integer, for example, 2 is 
defined as a basic number and supplied binary complex data 
is subjected to complex operation comprised of a single 
addition, subtraction and multiplication each so as to output 

35 binary complex data. 



Therefore, the N-point FFT is comprised of log 2 N-step 
stage and (N/2)log 2 N butterfly operations and an operation 
result of the FFT can be obtained depending on the number of 
multiplications of (2N)log 2 N times. Consequently, the 
5 N-point FFT is used as a Fourier transformation method ensuring 
a high operating efficiency. 

The FFT or discrete Fourier transformation (DFT) 
conducted in the above manner is usually used, so that data 
stream, sampled successively at a predetermined time interval 

10 and then supplied, is subjected to Fourier transformation. 
In that Fourier transformation, the sampled data stream is 
stored in a memory circuit in succession temporarily and when 
the number of data stored in the memory circuit temporarily 
reaches N, the Fourier transformation processing for that N 

15 data is started. 

During the Fourier transformation processing, data 
stored in the memory circuit temporarily is kept as it is and 
when the arithmetic operation is terminated, new data is 
supplied to the memory circuit again so that it is stored therein , 

20 Then, after the predetermined amount N of data is stored, the 
Fourier transformation processing is started again. 

However, because data is supplied without any 
interruption when the Fourier operation is executed also, a 
second memory circuit for accumulating N-data stream may be 

25 provided as another memory . Each time when N data are supplied, 
these memory circuits execute temporary storage processing 
and arithmetic operation processing alternately, so that the 
Fourier transformation is carried out on continuously supplied 
time-series data. 

30 In the Fourier transformation executed in such a manner , 

the arithmetic operation is carried out on N-data stream stored 
temporarily such that the data stream is treated as a block 
unit. Therefore, in such Fourier transformation processing, 
a delay of time corresponding to N-data stream is generated, 

35 so that a transformation result cannot be obtained at real 
time . 



Thus , although the data stream is supplied successively, 
the result of the Fourier analysis can be obtained only at 
N-samplings time interval. The delay time is not preferable 
for application of the Fourier analysis, thereby restricting 
5 its application field. 

In order to secure a result of the Fourier transformation 
each time when sampled data is supplied, by executing the 
Fourier transformation processing on N-point data containing 
new data supplied successively to be sampled, to solve the 

10 above -described problem, N-point Fourier transformation 
processings need to be carried out within a sampling period. 
Even in case where FFT operation method developed for such 
fast arithmetic operation processing is employed, N-time 
operating speed is demanded to secure a Fourier operation 

15 result computed continuously from the FFT. In a normal case, 
it is difficult to conduct such super high speed FFT to supply 
an obtained operation result . 

Although in the Fourier transformation, generally, the 
FFT operation processing method is used, the FFT operation 

20 processing is carried out on N-point data stream, which is 
sampled with a clock signal whose sampling frequency is fs, 
quant it ized and supplied, usually at a frequency interval of 
f s/N. 

In the arithmetic operation, data is fetched into other 
25 buffer memory even during the Fourier transformation 
processing, so that the Fourier analytic processing is carried 
out on time-series data supplied during the operation also. 
In a period in which data stored in one buff er memory temporarily 
is being processed, data supplied to the other memory is fetched 
30 thereinto. When fetching of N data is terminated, those 
arithmetic operation and data fetching processing are switched 
so as to execute the Fourier operation processing. However, 
this method needs two sets of the buffer memories and FFT 
operation processing means, which is not a favorable method 
35 in economic viewpoint . 

Further, because this method carries out block 



processing in which supplied N time- series data are treated 
as a block, the result of the Fourier transformation on the 
N-point data fetched in is not outputted until N sampling clock 
time passes, and the analysis result obtained at that time 
5 is only a result of the Fourier transformation processing for 
every N samples . 

The result of the Fourier transformation on the newest 
N-point data containing data, which is sampled successively 
and supplied, cannot be outputted at real time. If it is 

10 intended to output the transformation result at real time, 
the Fourier transformation processing in each sampling period 
is needed. However, if the Fourier transformation is carried 
out continuously at every sampling time interval, the amount 
of arithmetic operation per unit time becomes tremendous , which 

15 is not a realistic method. 

On the other hand, as a method for obtaining a result 
of continuous Fourier transformation processing, "Fourier 
transformation device and Fourier transformation method" has 
been proposed in Japanese Patent Application Laid-Open No. 

20 Hl-59454. 

This publication describes a method for 
Fourier- transforming a vibration waveform value sampled and 
supplied. According to the same method, a difference between 
a newly supplied vibration waveform value and an old vibration 

25 waveform value already supplied and used for Fourier 
transformation processing is obtained and a newer complex 
vibration value than the old complex vibration value provided 
by the arithmetic operation is obtained each time when a 
sampling value of the vibration waveform is supplied. 

30 Although frequency analysis may be considered as an 

application using the continuous Fourier transformation, if 
such application is considered, often any frequency band is 
analyzed at any resolution and then, that application is 
executed depending on the result . 

35 However, the aforementioned publication does not 

consider those points. 



Input data stream described in the same publication 
refers to only real number. Although if the data stream to 
be sampled is data stream which can be expressed in time series , 
it is permitted that the data stream is limited to the real 
5 number, imaginary number input signal as well as real number 
input signal may be sometimes used in Fourier operation for 
use in communication field, and therefore, supplied complex 
signals have to be subjected to transformation processing. 
However, the conventional fast Fourier transformation 
10 method, in which continuously supplied data is subjected to 
Fourier transformation to secure an operation result, is not 
capable of carrying out complex Fourier transformation on data 
stream supplied in the form of complex number. 

15 SUMMARY OF THE INVENTION 

Accordingly, an object of the present invention is to 
provide a continuous Fourier transformation method suitable 
for transforming any frequency band desired to be subjected 

20 to Fourier analysis at a desired resolution. 

Another object of the present invention is to provide 
a continuous Fourier transformation method enabling Fourier 
analysis on data stream supplied in the form of complex number 
and suitable for transforming any frequency band desired to 

25 be subjected to Fourier analysis at a desired resolution. 

To achieve the above object, there is provided a 
recursive discrete Fourier transformation method wherein data 

values x(t), x(t+l), x(t+2), x(t+3) x(t+N-l), x(t+N) 

sampled at times t, t + 1, t + 2 , t + 3 t+N-1, t+N (N is a 

30 positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which 
is degree k (kis 0 or a positive integer smaller than N) obtained 
by carrying out complex Fourier transformation on the data 

35 stream, is obtained such that a real part Xr(k, t) thereof 
and an imaginary part X ± (k, t) thereof are complex Fourier 
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coefficients, the method comprising: a first step of storing 
the data stream x(t) , x(t+l) , x(t+2) , x(t+3) , .... x(t+N-l) 
supplied since time t at time t+N-1 into a first memory means 
temporarily; a second step of obtaining the complex Fourier 
5 coefficients X r (k, t) and Xj.(k, t) of the data stream stored 
temporarily at the first step; a third step of storing the 
complex Fourier coefficients X r (k, t) and X ± (k, t) obtained 
at the second step into the second memory means temporarily; 
and a fourth step of by using data value x(t+N) supplied at 

10 time t+N, data value x(t) stored in the first memory means 
temporarily and the complex Fourier coefficients X r (k, t), 
Xi(k,t) stored in the second memory means temporarily, 
obtaining complex Fourier coefficients X r (k, t+1) and X ± (k, 
t+1) to the data stream supplied since time t+1 with respect 

15 to a positive constant value A for giving an amplitude value 
corresponding to a difference between the x(t+N) and the x(t) 
according to a following equation: 




20 

According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
transformation and then, the oldest sample value used for 
complex Fourier transformation is subtracted from the value 

25 of a new sample to be supplied next so as to obtain a subtraction 
value . By deleting the oldest sample at the same time , discrete 
Fourier operation result on new N samples can be obtained based 
on that subtraction value and the result of the preceding 
complex Fourier operation. Thus, different from the 

30 conventional Fourier operation that is carried out after N 
sample data is supplied, a Fourier operation result can be 



obtained quickly in a single sample period. 

Further, to achieve the above object, there is provided 
a recursive discrete Fourier transformation method wherein 
with sampling frequency as fs, data values x(t) , x(t+l) , x(t+2) » 
5 x(t+3), .... x(t+N-l), x(t+N) sampled at times t, t+1, t+2, 
t+3, .... t+N-1, t+N (N is a positive integer which is 1 or 
more) each having an equal interval are supplied and with such 
N data values supplied since time t as a data stream, a frequency 
interval given with a minimum frequency fl and a maximum 

10 frequency f2 to the data stream is regarded as a measuring 
frequency interval , a frequency interval obtained by dividing 
that measuring frequency interval by the N is assumed to be 
an analysis frequency interval and a result of the frequency 
analysis obtained by carrying out complex Fourier 

15 transformation at every analysis frequency interval is 
obtained as a frequency component which is k ( k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part X r (k, t) thereof and an imaginary part Xj_(k, 
t) thereof are complex Fourier coefficients, the method 

20 comprising: a first step of storing the data streamx( t ) , x(t+l) , 
x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; a second step 
of obtaining the complex Fourier coefficients X r (k, t) and 
Xi(k, t) of the data stream stored temporarily at the first 

25 step; a third step of storing the complex Fourier coefficients 
X r (k, t) and X ± (k, t) obtained at the second step into the 
second memory means temporarily; and a fourth step of by using 
data value x( t+N) supplied at time t+N, data value x( t ) stored 
in the first memory means temporarily and the complex Fourier 

30 coefficients X r (k, t), X ± (k,t) stored in the second memory 
means temporarily, obtaining complex Fourier coefficients 
X r (k, t+1) and Xi(k, t + 1) within the frequency interval with 
the minimum frequency f 1 and the maximum frequency f 2 to the 
data stream supplied since time t+1 with respect to a positive 

35 constant A for giving an amplitude corresponding to a 
difference between the x(t+N) and the x(t) according to 
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following equations: 



x, (K, + 1) - Jjr, (*.«)♦ -}At * n) - *(<)]}* cos j 2 | + /i] J 

(m + 1) - jt,M« cos jaf + fl]} - {*, M+^W + "M')| 



According to the present invention, data stream 
comprised of supplied N samples is subjected to complex Fourier 
transformation and then, the oldest sample value used for 

10 complex Fourier transformation is subtracted from the value 
of a new sample to be supplied next so as to obtain a subtraction 
value . By deleting the oldest sample at the same time , discrete 
Fourier operation result on newN samples can be obtained based 
on that subtraction value and the result of the preceding 

15 complex Fourier operation. At this time, a frequency band 
for analysis for the discrete Fourier operation is set up with 
f 1 and f2 preliminarily and the Fourier operation is carried 
out with a basic frequency whose resolution is set to N . Thus , 
different from a conventional case in which the frequency band 

20 and resolution are automatically set up when the Fourier 
operation is carried out after the N sample data is supplied, 
a result of the Fourier operation within any frequency band 
at any resolution can be obtained quickly in a single sampling 
period . 

25 According to a preferred embodiment of the present 

invention, inverse discrete Fourier transformation data 
subjected to inverse discrete Fourier transformation using 
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a positive constant a is supplied so as to obtain a complex 
Fourier transformation coefficient to the supplied inverse 
discrete Fourier transformation data, the discrete Fourier 
transformation being carried out using the positive constant 
5 A corresponding to the positive constant a. 

According to this embodiment, in addition to the above 
described effect, the operation can be carried out by FFT in 
which the positive constant A is set up by selecting 
appropriately 1, N or a square root of N. This is employed 

10 in, for example, communication apparatus. In case where 
information signal to be transmitted is subjected to inverse 
fast Fourier transformation (IFFT) processing to be converted 
to time-series signal and that time-series signal is 
transmitted, in Fourier operation method for obtaining 

15 information signal by supplying that transmitted signal to 
the FFT for executing the FFT operation, the constant value 
A corresponding to the constant value a, which is , for example , 
1, N or the square root of N used for the IFFT operation, is 
employed in the FFT operation processing circuit, so as to 

20 construct the FFT which operates in a complementary way with 
the IFFT to decode the information signal. By using the 
constant value A corresponding to the value a used for the 
IFFT operation for a signal supplied from other mating system, 
system having an excellent characteristic can be established. 

25 According to a preferred embodiment of the present 

invention, k indicating a degree of the complex Fourier 
coefficient is set to a desired value and Fourier 
transformation to the set value k is repeated at every time 
having an equal interval. 

30 According to this embodiment, in addition to the above 

described effect, particularly by setting up the value of k 
to a specific value, a result of the Fourier analysis on a 
specific frequency determined with the value k can be computed 
quickly. Thus, if there is a system controlled according to 

35 the result of the Fourier transformation on that specific 
frequency, the control of that system can be carried out at 
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a minimum delay time, so that response characteristic of the 
entire system can be improved. Thus, system having a better 
quality can be constructed. 

Further, to achieve the above object, there is provided 
5 a recursive discrete Fourier transformation method wherein 

datavaluesx(t) , x(t+l) , x(t+2) , x(t+3) x(t+N-l) , x(t+N) 

sampled at times t, t+1, t+2, t+3, t+N-1, t+N (N is a 

positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 

10 since time t as a data stream, a frequency component, which 
is degree k (kis 0 or a positive integer smaller thanN) obtained 
by carrying out complex Fourier transformation on the data 
stream, is obtained such that a real part X r (k, t) thereof 
and an imaginary part X±(k, t) thereof are complex Fourier 

15 coefficients, the method comprising: a first step of storing 
the data stream x(t) , x(t+l) , x(t+2) , x(t+3) , ...,x(t+N-l) 
supplied since time t at time t+N-1 into a first memory means 
temporarily; a second step of obtaining the complex Fourier 
coefficients X r (k, t) and Xi(k, t) of the data stream stored 

20 temporarily at the first step; a third step of storing the 
complex Fourier coefficients X r (k, t) and X ± (k, t) obtained 
at the second step into the second memory means temporarily; 
and a fourth step of by using the complex Fourier coefficients 
X r (k, t) - jXi(k, t) stored in the secondmemory means temporarily , 

25 obtaining complex Fourier coefficients X r (k, t+1) - jXi(k, 
t+1) to the data stream supplied since time t+1 based on a 
transfer function expressed in a following equation. 



H (z) = a{i-z- n \ 



. . \ n Jik 1 
cos 2 — - ; sin 2 — 
# \ N\ 



l-2cos 2 — z" 1 +z" 



where A is a positive constant for providing [x(t+N) - x(t)] 
30 with an amplitude. 

According to the present invention, data stream 
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comprised of supplied N samples is subjected to complex Fourier 
transformation and then, the oldest sample value used for 
complex Fourier transformation is subtracted from the value 
of a new sample to be supplied next so as to obtain a subtraction 
5 value. By deleting the oldest sample at the same time, discrete 
Fourier coefficients to new N samples can be obtained by 
operation with a predetermined transfer function based on that 
subtraction value and the preceding complex Fourier operation 
result. Thus, different from the conventional Fourier 

10 operation that is carried out after N sample data is supplied, 
a Fourier operation result can be obtained quickly in a single 
sample period. 

Further, to achieve the above object, there is provided 
a recursive discrete Fourier transformation method wherein 

15 with sampling frequency as f s , data values x(t), x(t+l),x(t+2), 
x(t+3), x(t+N-l), x(t+N) sampled at times t, t+1, t+2, 

t+3, . . . , t+N-1, t+N (N is a positive integer which is 1 or 
more) each having an equal interval are supplied and with such 
N data values supplied since time t as a data stream , a frequency 

20 interval given with a minimum frequency fl and a maximum 
frequency f2 to that data stream is regarded as a measuring 
frequency interval , a frequency interval obtained by dividing 
that measuring frequency interval by the N is assumed to be 
an analysis frequency interval and a result of the frequency 

25 analysis provided by carrying out complex Fourier 
transformation at every analysis frequency interval is 
obtained as a frequency component which is k (k is 0 or a positive 
integer smaller than N) times the frequency interval, such 
that a real part X r (k, t) thereof and an imaginary part Xi(k, 

30 t) thereof are complex Fourier coefficients, the method 
comprising: a first step of storing the datastreamx(t) , x(t+l) , 
x(t+2), x(t+3), x(t+N-l) supplied since time t at time 

t+N-1 into a first memory means temporarily; a second step 
of obtaining the complex Fourier coefficients X r (k, t) and 

35 Xi(k, t) of the data stream stored temporarily at the first 
step; a third step of storing the complex Fourier coefficients 
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X r (k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and a fourth step of by using 
the complex Fourier coefficient X r (k, t) - jXi(k, t) stored 
in the second memory means temporarily, obtaining complex 
5 Fourier coefficients X r (k, t+1) - jXj. (k, t+1) at the frequency 
interval given by the minimum frequency fl and the maximum 
frequency f2 to the data stream supplied since time t+1 based 
on a transfer function expressed in a following equation: 

H(z) = A(l -z' n )j cosM-;sin[2^]-z-M 
w v y | l-2cos[2^?>- 1 +z~ 2 J 

10 where A is a positive constant for providing [x(t+N) - x(t)] 
with an amplitude value, and 

According to the present invention, data stream 

15 comprised of supplied N samples is sub jected to complex Fourier 
transformation and then, the oldest sample value used for 
complex Fourier transformation is subtracted from the value 
of a new sample to be supplied next so as to obtain a subtraction 
value . By deleting the oldest sample at the same time , discrete 

20 Fourier operation result on new N samples can be obtained based 
on that subtraction value and the result of preceding the 
complex Fourier operation. At this time, a frequency band 
for analysis for the discrete Fourier operation is set up with 
f 1 and f 2 preliminarily and that Fourier operation is carried 

25 out with a basic frequency whose resolution is set to N . Thus, 
different from a conventional case in which the frequency band 
and resolution are automatically set up when the Fourier 
operation is carried out after the N sample data is supplied, 
a result of the Fourier operation within any frequency band 

30 at any resolution can be obtained quickly in a single sampling 
period. 

According to a preferred embodiment of the present 
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invention, the positive constant value A for providing with 
an amplitude corresponding to a difference between the x(t+N) 
and the x( t ) is capable of being set selectively with an inverse 
number of square root of 1, N or 1/N. 
5 According to this embodiment, in addition to the above 

described effect , the operation can be carried out by FFT in 
which the positive constant A is set up by selecting 
appropriately 1, 1/N or an inverse number of a square root 
of N. This is employed in, for example, communication 

10 apparatus . In case where information signal to be transmitted 
is subjected to inverse fast Fourier transformation (IFFT) 
processing to be converted to time- series signal and that 
time- series signal is transmitted, in Fourier operation method 
for obtaining information signal by supplying that transmitted 

15 signal to the FFT for executing the FFT operation, the constant 
value A corresponding to the constant value a, which is, for 
example, 1, 1/N or the inverse number of the square root of 
N used for the IFFT operation, is employed in the FFT operation 
processing circuit, so as to construct the FFT which operates 

20 in a complementary way with the IFFT to decode the information 
signal. By using the constant value A corresponding to the 
value a used for the IFFT operation for a signal supplied from 
other mating system, systemhaving an excellent characteristic 
can be established. 

25 According to a preferred embodiment of the present 

invention, k indicating the degree of the complex Fourier 
coefficient is set to a desired value and the Fourier 
transformation to the set value k is repeated at every time 
having an equal interval . 

30 According to this embodiment , in addition to the above 

described effect, particularly by setting up the value of k 
to a specific value, a result of the Fourier analysis on a 
specific frequency determined with the value k can be computed 
in a short time. Thus, if there is a system controlled 

35 according to the result of the Fourier transformation on that 
specific frequency, the control of that system can be carried 
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out at a minimum delay time, so that response characteristic 
of the entire system can be improved. Thus, system having 
a better quality can be constructed. 

Further, to achieve the above object, there is provided 
5 a recursive discrete Fourier transformation method for 
obtaining complex Fourier coefficients, wherein data values 
x r (t) + jXi(t) , x x (t+l) + jXiCt + l), x r (t+2) + jXi(t+2), x r (t + 3) 
+ jXi(t+3), . . . ,-x r (t+N-l) + 3Xi(t+N-l), x r (t+N) + jXi(t+N) 
sampled at times t, t+1, t+2, t+3, .... t+N-1, t+N (N is a 

10 positive integer which is 1 or more) each having an equal 
interval are supplied and with such N data values supplied 
since time t as a data stream, a frequency component, which 
is degree k ( k is 0 or a positive integer smaller than N ) obtained 
by carrying out complex Fourier transformation on the data 

15 stream, is obtained such that a real part X r (k, t) thereof 
and an imaginary part Xi(k, t) thereof are complex Fourier 
coefficients, the method comprising: a first step of storing 
the data stream x r ( t )+jXi(t) , x r (t+l) + jXi(t+l), x r (t+2) + 
3Xi(t+2), x r (t+3) + 3Xi(t+3), x r ( t+N-1) + jXi(t+N-l) 

20 supplied since time t at time t+N-1 into a first memory means 
temporarily; a second step of obtaining the complex Fourier 
coefficients X r (k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; a third step of storing the 
complex Fourier coefficients X r (k, t) and Xi(k, t) obtained 

25 at the second step into the second memory means temporarily; 
and a fourth step of by using the data value x r (t+N) + jXi(t+N) 
supplied at time t+N, the data value x r (t) + jxi(t) stored 
in the first memory means temporarily and the complex Fourier 
coefficient X r (k, t) and Xi(k, t) stored in the second memory 

30 means temporarily and used recursively, obtaining complex 
Fourier coefficients X r (k, t+1) and X ± (k, t + 1) to the data 
stream supplied since time t+1 with respect to a positive 
constant A for giving an amplitude value corresponding to the 
difference between the x r (t+N) and the x r (t) according to a 

35 following equation: 
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X. t {k,t + 1) = jz . (it, 0 + -jfo (t + N) - x, (* )]j cos J^J 
+ {2T,M + ^kfr+^)-^W]}«in[2f] 

5 

According to the present invention, data stream 
comprised of supplied N complex data is subjected to Fourier 
transformation, so that a Fourier operation result to new N 
complex data can be obtained based on the preceding Fourier 

10 transformation result, data value of new complex data to be 
supplied next and the oldest complex data used for executing 
the Fourier transformation. Therefore, different from the 
conventional case where the Fourier operation is carried out 
after the N sample data is supplied, the Fourier operation 

15 result can be obtained quickly in the period of a single complex 
data. 

Further, to achieve the above object, there is provided 
a recursive discrete Fourier transformation method for 
obtaining complex Fourier coefficient wherein with sampling 

20 frequency as fs, data values x r (t)+jXi(t) , x r (t+l) + jXi(t+l), 
x r (t+2) + jx ± (t+2), x r (t+3) + jXi(t+3), x r (t+N-l) + 

ZjXiCt+N-l), x r (t+N) + jXi(t+N) sampled at times t, t+1, t + 2, 
t+3, . . . , t+N-1, t+N (N is a positive integer which is 1 or 
more) each having an equal interval are supplied and with such 

25 N complex data values supplied since time t as a data stream, 
a frequency interval given with a minimum frequency fl and 
a maximum frequency f 2 to that data stream is regarded as a 
specified frequency interval, a frequency interval obtained 
by dividing the specif led frequency interval by the N is assumed 
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to be a minimum frequency interval and a result of the Fourier 
transformation provided by carrying out complex Fourier 
transformation at every minimum frequency interval is obtained 
as a frequency component which is k (k is 0 or a positive integer 
5 smaller than N) times the minimum frequency interval, such 
that a real part X r (k, t) thereof and an imaginary part X ± (k, 
t) thereof are complex Fourier coefficients, the method 
comprising: a first step of storing the data stream x r (t) + 
jXi(t), x r (t+l) + jXi(t+l), x x (t+2) + jx ± (t+2), x x (t+3) + 

10 jx ± (t+3), x r (t+N-l) + jxi(t+N-l), x r (t+N) + jXi(t+N) 

supplied since time t at time t+N-1 into a first memory means 
temporarily; a second step of obtaining the complex Fourier 
coefficients X r (k, t) and Xi(k, t) of the data stream stored 
temporarily at the first step; a third step of storing the 

15 complex Fourier coefficients X r (k, t) and X ± (k, t) obtained 
at the second step into the second memory means temporarily; 
and a fourth step of by using the data value x r (t+N) + jXi(t+N) 
supplied at time t+N, the data value x r (t) + jx ± (t) stored 
temporarily in the first memory means and complex Fourier 

20 coefficients X r (k, t) and Xi(k, t) stored temporarily in the 
second memory means and used recursively, obtaining complex 
Fourier coefficient X r (k, t+1) and Xi(k, t+1) in the specified 
frequency interval given by the minimum frequency f 1 and the 
maximum frequency f2 to the data stream supplied since time 

25 t + 1 withrespect to apositive constant Af or giving an amplitude 
value corresponding to a difference between the x(t+N) and 
the x(t) according to a following equation: 




(f2-n)k 



N-l 




(f2-fl)k 



1) 



N-l 



30 




• (/2-/1)* 



+ /1 



N-l 
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According to the present invention, data stream 
comprised of supplied N complex data is subjected to Fourier 
5 transformation and then , when it is intended to obtain a result 
of the Fourier operation to new N complex data on the basis 
of the result of the preceding Fourier transformation, data 
value of new complex data to be supplied next and the oldest 
complex data used for Fourier, the frequency band to be 

10 Fourier -transformed is set up with fl and f2 preliminarily 
and further, the Fourier operation is carried out using the 
basic frequency set with the resolution N. Therefore, 
different from the conventional case where the frequency band 
and the resolution are automatically determined when the N 

15 sample data are supplied for the Fourier operation, a Fourier 
operation result to any frequency band and resolution can be 
obtained quickly and easily in the period of a single sample. 

According to a preferred embodiment of the present 
invention, by using data value y r (t ) + jyi(t) supplied at time 

20 t, data value y r (t+N) + jyi(t+N) supplied at time t+N and complex 
inverse discrete Fourier coefficients Y r (k, t) and Yi(k, t) 
with the real part and the imaginary part obtained with respect 
to Ndata values y r (t) + 3y±(t) , y r ( t+1) + jyi(t+l) , . . . , y r (t+N-l) 
+ jyx(t+N-l) supplied from time t to time t+N-1, a real part 

25 Y r (k, t+1) and an imaginary part Yi(k, t+1) of each complex 
inverse discrete Fourier coefficient of N data values supplied 
since time t+1 are obtained with respect to a positive constant 
value B for giving an amplitude corresponding to a difference 
between the y r (t+N) and the y r (t) as inverse discrete Fourier 

30 transformation data according to following equations, 

y r (K, + 1) - {y, (*, 0 + |[y, (< + N) - y, (r )]J cos [ 2 £] 
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♦ {r,M+£U(<+*)-»M]}™[2£] 

and after the obtained inverse discrete Fourier transformation 
5 data is supplied, discrete Fourier transformation on the 
supplied inverse discrete Fourier transformation data is 
carried out, the recursive discrete Fourier transformation 
conducting discrete Fourier transformation by using the 
positive constant value A corresponding to the positive 

10 constant value B. 

According to this embodiment , in addition to the above 
described effect, the Fourier transformation can be carried 
out such that the positive constant A is set up to 1, N or 
a square root of N. This is employed in, for example, 

15 communication apparatus. In case where information signal 
to be transmitted is subjected to inverse fast Fourier 
transformation (IFFT) processing to be converted to 
time-series signal and that time-series signal is transmitted, 
in the Fourier operation method for obtaining information 

20 signal by supplying that transmitted signal to the FFT for 
executing the FFT operation , the constant value A corresponding 
to the constant value B, which is, for example, 1, N or the 
square root of N used for the IFFT operation, is employed in 
the FFT operation processing circuit, so as to construct the 

25 FFT which operates in a complementary way with the IFFT. By 
using the constant value A corresponding to the value B used 
for the IFFT operation for a signal supplied from other mating 
system, system containing Fourier transformation processing 
having an excellent characteristic can be established. 

30 According to a preferred embodiment of the present 

invention, k indicating the degree of the complex Fourier 
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coefficient is set to a desired value and the Fourier 
transformation to the set value k is repeated at every time 
having an equal interval. 

According to this embodiment, in addition to the above 
5 described effect, particularly by setting up the value of k 
to a specific value, the Fourier analysis on a specific 
frequency determined with the value k can be carried out in 
a short time. Thus, if there is a system controlled according 
to the result of the Fourier transformation on that specific 
10 frequency, the control of that system can be carried out at 
a minimum delay time, so that response characteristic of the 
entire system can be improved. Thus, system whose response 
characteristic is a short time, having a better quality, can 
be constructed. 

15 Further, to achieve the above object, there is provided 

a recursive inverse discrete Fourier transformation method 
for obtaining complex inverse Fourier coefficient wherein 
complex data values x r (t) +3Xi(t), x r (t+l) + jXi(t+l), x r (t+2) 
+ JXi(t+2), x r (t+3) + jx ± (t+3), x r (t+N-l) + jXi(t+N-l), 

20 x r (t+N) + jXi(t+N) sampled at times t, t+1, t+2, t+3 

t+N-l, t+N (N is a positive integer which is 1 or more) each 
having an equal interval are supplied and with such N data 
values supplied since time t as a data stream , an inverse complex 
Fourier transformation component, which is degree k (k is 0 

25 or a positive integer smaller than N) obtained by carrying 
out inverse complex Fourier transformation on the data stream, 
is obtained such that a real part X r (k, t) thereof and an 
imaginary part Xi(k, t) thereof are inverse complex Fourier 
coefficients, the method comprising: a first step of storing 

30 the data stream x r (t) + jXi(t) , x r (t+l) + jxi(t + l), x r (t+2) + 
jXi(t+2), x r (t+3) + jXi(t+3), x r (t+N-l) + jx±(t+N-l) 

supplied since time t at time t+N-l into a first memory means 
temporarily; a second step of obtaining the inverse complex 
Fourier coefficients X r (k, t) and Xi(k, t) of the data stream 

35 stored temporarily at the first step; a third step of storing 
the inverse complex Fourier coefficients X r (k, t) and Xi(k, 
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t) obtained at the second step into the second memory means 
temporarily ; and a fourth s tep of by using the data value x r ( t +N ) 
+ jXi(t+N) supplied at time t+N, the data value x r (t) + jXi(t) 
stored in the first memory means temporarily and the inverse 
5 complex Fourier coefficient X r (k, t) and Xi(k, t) stored in 
the second memory means temporarily, obtaining inverse complex 
Fourier coefficients X r (k, t+1) and Xi(k, t + 1) to the data 
stream supplied since time t+1 with respect to a positive 
constant B for giving an amplitude value corresponding to the 
10 difference between the x r (t+N) and the x r (t) according to a 
following equation: 



According to the present invention, data stream 
comprised of supplied N complex data is subjected to inverse 
Fourier transformation, so that an inverse Fourier operation 
result to new N complex data can be obtained based on the 

20 preceding inverse Fourier transformation result, data value 
of new complex data to be supplied next and the oldest complex 
data used for executing the inverse Fourier transformation. 
Therefore, different from the conventional case where the 
inverse Fourier operation is carried out after the N sample 

25 data is supplied, the inverse Fourier operation result can 
be obtained quickly in the period of a single complex data. 






15 




Further, to achieve the above object, there is provided 
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a recursive inverse discrete Fourier transformation method 
for obtaining complex inverse Fourier coefficients wherein 
with sampling frequency as fs, data values x r (t) + jXi(t), 

X r (t+1) + jXi(t + l) , X r (t + 2) + jXi(t + 2) . x r (t + 3) + jXi(t + 3) , 

5 x r (t+N-l) + jXi(t+N-l), x r (t+N) + jXi(t+N) sampled at times 
t, t+1, t+2, t+3, t+N-1, t+N (N is a positive integer 

which is 1 or more) each having an equal interval are supplied 
and with such N data values supplied since time t as a data 
stream, a frequency interval given with a minimum frequency 

10 f 1 and a maximum frequency f 2 to the data stream is regarded 
as a specified frequency interval, a frequency interval 
obtained by dividing that specified frequency interval by the 
N is assumed to be a minimum frequency interval and a result 
of the inverse Fourier transformation provided by carrying 

15 out inverse complex Fourier transformation at every minimum 
frequency interval is obtained as a frequency component which 
is k (k is 0 or a positive integer smaller than N) times the 
minimum frequency interval, in the form of a real part X r (k, 
t) and an imaginary part Xi(k, t), the method comprising: a 

20 first step of storing the data stream x r (t)+jXi( t) , x r (t+l) 
+ jXi(t+l), x r (t+2) + jXi(t+2), x r (t+3) + 3Xi(t+3), 
x r (t+N-l) + jXi(t+N-l), x r (t+N) + jx±( t+N) supplied since time 
t at time t+N-1 into a first memory means temporarily; a second 
step of obtaining the inverse Fourier coefficients X r (k, t) 

25 and Xi(k, t ) of the data stream stored temporarily at the first 
step; a third step of storing the inverse Fourier coefficients 
X r (k, t) and Xi(k, t) obtained at the second step into the 
second memory means temporarily; and a fourth step of by using 
the data value x r (t+N) + jXi(t+N) supplied at time t+N, the 

30 data value x r (t) + jxi(t) stored temporarily in the first memory 
means and preceding inverse Fourier coefficients X r (k, t) and 
Xi(k, t) stored temporarily in the second memory means, 
obtaining inverse Fourier coefficient X r (k, t+1) and Xi(k, 
t+1) in the specified frequency interval given by the minimum 

35 frequency f 1 and the maximum frequency f 2 to the data stream 
supplied since time t+1 with respect to a positive constant 
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B for giving an amplitude value corresponding to a difference 
between the x(t+N) and the x(t) according to following 
equations : 

X, (k,, + 1) - {*, ±k (' + *)-*, Mljxcos )2^[%^ + /l] J 
+ {z,(M) + I[,(,^)-,X0]}«in{2|[%i^ + /l]} 

X,(Kt + 1) - {jr,M + }tx,fr + *)-*,W]}xa»[2i[i£^i + a]} 
-jx,fc<) + |k(< + ^)-^W])xs i n{2i[I^)* + / 1 ]J 



According to the present invention, data stream 

10 comprised of supplied N complex data is subjected to inverse 
Fourier transformation and then, when it is intended to obtain 
a result of the inverse Fourier operation to new N complex 
data on the basis of the result of the preceding inverse Fourier 
transformation, data value of new complex data to be supplied 

15 next and the oldest complex data used for inverse Fourier, 
the frequency band to be inverse-Fourier-transformed is set 
up with f 1 and f 2 preliminarily and further , the inverse Fourier 
operation is carried out using the basic frequency set with 
the resolution N . Therefore , different from the conventional 

20 case where the frequency band and the resolution are 
automatically determined when the N sample data are supplied 
for the inverse Fourier operation , an inverse Fourier operation 
result to any frequency band and resolution can be obtained 
quickly and easily in the period of a single sample. 

25 According to a preferred embodiment of the present 

invention, k indicating the degree of the complex Fourier 
coefficient is set to a desired value and the inverse Fourier 
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transformation to the set value k is repeated at every time 
having an equal interval. 

According to this embodiment , in addition to the above 
described effect, particularly by setting up the value of k 
5 to a specific value, the inverse Fourier analysis on a specific 
frequency determined with the value k can be carried out in 
a short time. Thus, if there is a system for generating a 
signal according to the result of the inverse Fourier 
transformation on that specific frequency, an operation of 

10 that system can be executed at a minimum delay time, so that 
response characteristic of the entire system can be improved. 
Thus , because the response characteristic of the entire system 
can be improved, system whose response characteristic is 
shorter, having an excellent quality, can be constructed. 

15 The nature, principle and utility of the invention will 

become more apparent from the following detailed description 
when read in conjunction with the accompanying drawings . 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

In the accompanying drawings : 

Fig.l is a diagram showing schematically a recursive 
discrete Fourier transformation device according to an 
embodiment of the present invention; 

25 Fig. 2 is a diagram showing schematically a relation 

between data value obtained by sampling the waveform of signal 
supplied to the recursive discrete Fourier transformation 
device of the embodiment of the present invention and a 
corresponding DFT operation; 

30 Fig . 3 is a diagram showing the structure of the recursive 

discrete Fourier transformation , which carries out arithmetic 
operation using a transfer function according to the embodiment 
of the present invention; 

Fig. 4 is a diagram showing an analysis frequency band 

35 in the recursive discrete Fourier transformation device of 
the embodiment of the present invention and a relation of the 
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frequency relating to setting of resolution; 

Fig. 5 is a diagram showing the structure for an analysis 
frequency band and the setting of the resolution in the 
recursive discrete Fourier transformation device which 
5 carries out arithmetic operation using a transfer function 
of the embodiment of the present invention; and 

Fig. 6 is a diagram showing a schematic structure of the 
recursive inverse discrete Fourier transformation device of 
the embodiment of the present invention. 

10 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 



Hereinafter, the preferred embodiment of the recursive 
discrete Fourier transformation method and the recursive 
15 inverse discrete Fourier transformation method of the present 
invention will be described with reference to the accompanying 
drawings . 

Fig . 1 shows an embodiment of a recursive discrete Fourier 
transformation device applying the recursive discrete Fourier 

20 transformation method, which will be described below with 
reference to the same Figure. 

The recursive discrete Fourier transformation device 
comprises a data updating portion 1 to which data sampled in 
a predetermined time interval is supplied and which stores 

25 temporarily supplied new N (N: positive integer ) data therein, 
a basic frequency setting portion 2 in which the basic frequency 
for carrying out the discrete Fourier transformation is set 
up, a recursive DFT operating portion 3 for carrying out 
recursive DFT operation, and a memory 4 for storing temporarily 

30 computed data. 

Next , an operation of the recursive discrete Fourier 
transformation device having such a structure will be 
described. 

First, supplied data is sampled with a sampling circuit 
35 (not shown) every predetermined time interval and then, the 
sampled and quantitized discrete data is supplied to the data 
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updating portion 1 . 

The sampling circuit samples data values supplied at 
every time of a predetermined interval, t, t+1, t+2, t+3, 
t+N-1, t+N (N: natural number ) and generates data streams x( t ) , 
5 x(t+l), x(t+2), x(t+3), x(t+N-l), x(t+N) by converting 

the data values supplied at the respective times to sampling 
values corresponding to the respective times . 

The operation of the sampling circuit is the same as 
that of an A/D converter which converts supplied analog signal 
10 to digital signal, converting an analog signal voltage value 
supplied at a predetermined time interval, given in the form 
of an inverse number of a sampling frequency, to a digital 
signal value . The converted digital signal value is a voltage 
having an analogous relation with an analog voltage, which 
15 gives a pulse amplitude modulation signal or that voltage value 
is expressed with a binary digital value. 

In this manner, sampling data x(t) sampled at time t, 
x(t + l) sampled at time t + 1, ... ,-x( t+N-1) sampled at t+N-1, 
x(t+N) sampled at time t+N, . . . , are supplied to the data 
20 updating portion 1. 

The data-updating portion 1 updates the newest N data 
(N: positive integer) of the supplied data, and stores them 
temporarily. 

That is, data supply is started from x(t) and when data 
25 x(t+l), x(t+2), x(t+3) are supplied, all the supplied data 
x(t), x(t+l), x(t+2), x(t+3) are stored temporarily. This 
temporary storage operation is continued until data x( t+N-1) 
is entered and when data x( t+N-1) is entered, the total data 
number becomes N, so that data region of the data updating 
30 portion 1 becomes full. 

When next data x ( t +N ) is supplied with such a condition , 
the total data number becomes N+l . Thus, the data updating 
portion 1 subtracts x(t) from x(t+N) and then supplies data 
obtained by that subtraction to the recursive DFT operating 
35 portion 3 while removing the oldest data x(t) from the memory 
portion. Consequently, the data updating portion 1 
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temporarily stores Ndatax(t+1 ) , x(t+2) , x(t+3) ..... x(t+N-l) . 
x(t+N) . 

Likewise, when next data x(t+N+l) is supplied, the data 
updating portion 1 obtains x( t+N+l)-x(t+l) and supplies it 
5 to the recursive DFT operating portion 3 while removing x( t+1 ) 
from the memory. Consequently, the newest N data of supplied 
data is always stored in the data-updating portion 1 
temporarily. 

The N data stored temporarily is supplied to the 
10 recursive DFT operating portion 3 and a last FFT operation 
result stored in the memory 4 temporarily is supplied to the 
recursive DFT operating portion 3 as recursive data. Then, 
the recursive DFT operating portion 3 carries out recursive 
discrete Fourier transformation according to frequency 
15 resolution information set up by the basic frequency setting 
portion 2 in a manner described later and outputs a result 
of its arithmetic operation. 

Next, this recursive discrete Fourier transformation 
operation method will de described in detail with reference 
20 to a conventional Fourier operation as well. 

Fig . 2 explains schematically a relation between data 
value obtained by sampling the waveform of supplied signal 
in a sampling period ts and corresponding DFT operation. 

The same Figure shows a set of N real number data values 
25 x(t), x(t+l), x(t+2), x(t+3), x(t+N-l) sampled in a 

predetermined sampling period from time t and a set of N real 
number data values x(t+l), x(t+2), x(t+3), x(t+N-l), 
x(t+N) sampled from time t+1. 

Then, the value X(k, t) of discrete Fourier 
30 transformation obtained with respect to data stream which is 
N real number data values sampled since time t is defined in 
the following expression. 

X(k,t) = ^24nhxp^-j2(n-t)^ - (2) 



^|4 n )cos[2(„- i )f]-/^'fWs i n[2( n - ! )f] 

= X r (k,t)-jX i {k,t) 

, for k = 0, 1, N-l 

5 Transformation of real part X r (k, t) and imaginary part 

X ± (k, t) of complex Fourier transformation of sampling data 
stream x(t) ~ x(t+N~l) supplied since arbitrary time t is 
defined as follows. 



X,.{k,t) = 4= ,+ y x (n)cos\2{n-t)— 1 - (3) 

-IN [ N\ 




for k = 0, 1 N-l 



The transformation equation for the sampling data 
supplied since time t is defined in the above manner. Next, 

15 a transformation equation for a sampling data stream to be 
supplied since time t+1 will be described below. 

That is, the sampling data stream supplied since the 
time t+1 is expressed as a set of x(t + l) , x(t+2), x(t+3), 
•x(t+N-l), x(t+N). Because data x(t) is deleted from data 

20 stream supplied since time t while new data is added, 
transformation of the real part is expressed in a following 
equation and that expressed equation is developed as follows. 
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where , 

r r (^ + l)-^2*(»)oos|2(i»-/)^J - (6) 

Yfat +l) = -^2 x ( ?z ) sin [ 2 ("- / )y] "' (7) 

Then, if the oldest data x(t) is deleted and the newest 
5 datax(t+N) is incorporated, Y r (k, t+1) can be expressed with 
X r (k, t) and likewise, Yi(k, t+1) can be also expressed with 
Xi(k, t). 

As a result, the equation (6) and equation (7) are 
expressed as follows. 

10 Y,(k,t+l) = X r {k,t) + -^=[x(t + N)-x(t)] (8) 

y;(*,t + i) = jjt i .(jm) — (9) 



Next , the transformation equation for the imaginary part 
is also developed in the same way. 

That is, when new data x(t+N) is supplied, the 
15 transformation of the imaginary part is also developed as shown 
in a following equation. 

zA^+l) = ^|x(n)sin|2(/2-/-l)^j - (10) 

-Y t {k,t + l)cos [2^J -Y r (k,t + l)sin \lf j 

20 Therefore, the equations (5), (10) can be expressed based 

on the relation of the equations (8), (9). 

X r (k,t + 1) = \x r (k,t) + -~[x(t + N)-x(t)]x cos |2^J + X, (k,t)sin | 2 ^] J < 11 > 
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X i {k,t + l) = X i (k,t)cos ^Y^X r (k,t) + ^[x(t + N)-x(t)^sm j^J (12) 

In this way, the discrete Fourier transformation at time 
t+1 can be executed according to the equations (11), (12) by 
using the result of the discrete Fourier arithmetic operation 
5 at time t recursively and further using data value which is 
a difference between newly fetched sample value x(t+N) and 
a sample value x(t) to be deleted. 

It means that a result of Fourier transformation on the 
newest N-point data containing new data obtained by successive 
10 sampling at the sampling time interval shown in Fig. 2 can be 
obtained. 

Here, if it is regarded that the Fourier transformation 
executed in the above- describedmanner is a linearly invariable 
system in which continuous sampling data is supplied to a system 

15 having a certain coefficient and converted, this conversion 
processing can be regarded to be filtering processing. 

Then, if l/( square root of N) which gives an amplitude 
is assumed to be A and the real part X r (k,t) and the imaginary 
part Xi ( k , t ) of Fourier coefficient with respect to data stream 

20 supplied since time t are assumed to be y r , k (t) and yi, k (t) 
respectively, the aforementioned equations (11), (12) are 
treated as a transformation equation of output y to input x, 
so that the transformation equation of the Fourier coefficient 
of data supplied since time t+1 can be expressed as follows. 

25 y nk (t + 1) = {y r , (t) + A[x(t + N)-x(t )]}x cos + y iJc (t) sin | 2 ^J -(13) 

>U^ = >U0cos^j4^ -(14) 

(where A = ) 
Further, if these equations (13), (14) are subjected 
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to Z transformation, a result of the transformation can be 
expressed as shown in following equations (15), (16). 

Y r * M = d)*- 1 +A[x(z)z» - *(z)}x cos |2^J + Y Uk (z^ 1 sin [2^] -(15) 




-(16) 



5 If a transfer function of output Y(z) to input X(z) is 

assumed to be H(z), H(z) is given in a following equation. 

H(z) = Y(z)/X(z) 

Then, if the equations (15), (16) are treated as 
simultaneous equations and it is intended to obtain each 
10 transfer function of the real part and imaginary part, they 
can be obtained as shown in following equations (17) , (18) . 



H r (z) = A(l-z- N ). 



|l-2cos|2^|z- 1 + z- 2 



- (17) 



H i (z) = A(l-z' N ). 



-sin 2- 



l-2cos|2^jz- 1 +z- 2 



Like this, a transfer function indicating the recursive 
15 DFT is expressed in a following equation (19). 



H(z) = A[l~z- N 



cos 2 — - j sin 2 — 
1 ' N \ J N 



l-2cos[2— ]z- 5 +z- 2 
N 



(19) 
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In this equation, 1-Z" N indicates deletion of the oldest 
data upon updating the input data and a portion expressed with 
a fraction is so constructed that its numerator and its 
denominator are IIR filters (cyclic filter) having primary 
5 and secondary degrees respectively. 

Although the DFT is a FIR filter for carrying out 
processing to a finite point number N, if it is achieved with 
a IIR filter as shown in the aforementioned equation, it can 
be constructed with a simpler hardware than the FIR filter. 
10 Fig. 3 shows the structure for carrying out the discrete 

Fourier transformation achieved with the IIR filter. 

In the same Figure, data sampled at a predetermined 
sampling clock, for example, data x(t+N) sampled at time t+N 
is supplied to a Fourier operating portion 5. Here, the 
15 supplied sampling data is computed by the characteristic of 
a transfer function indicated asH(z), so that X r (k, t+1) and 
jX ± (k, t+1) are supplied as real part output and imaginary 
part output of the Fourier coefficient. 

The Fourier operating portion 5 having the transfer 
20 function H(z) indicated by the equation (19) carries out IIR 
filtering operation in succession for every supplied sampling 
data and supplies results of the operation of the DFT 
transformation corresponding to the real part and the imaginary 
part . 

25 That is, if data sampled at a certain sampling time t + 1 

is supplied to the Fourier operating portion 5 , complex Fourier 
coefficient X r (k. t + 1) - jX ± (k, t+1) corresponding to N sampling 
data supplied since the sampling time t+1 can be obtained. 

As described above, the Fourier operating portion 5 can 

30 carry out Fourier transformation processing on continuously 
supplied data stream continuously at the sampling time interval. 
Then, a result of arithmetic operation can be obtained by 
changing analysis resolution in the Fourier transformation, 
so that high-resolution frequency analysis can be carried out 

35 in time direction by changing the length N of data stream 
supplied, for example, increasing the value of N. 
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As for the change of the analysis resolution, in case 
of adaptive signal processing to be carried out using a result 
of analysis by, for example, the DFT operation or adaptive 
control for control system, which is an object, the discrete 
5 Fourier transformation needs to be analyzed to arbitrary 
frequency band at an arbitrary resolution (frequency 
interval) . 

Fig. 4 shows the relation of frequency, which should be 
obtained by an object discrete Fourier analysis. 

10 The same Figure shows recursive discrete Fourier 

transformation in which the minimum frequency of a frequency 
band to be analyzed is f 1 while the maximum frequency thereof 
is f2 and N-point resolution is provided between fl and f 2 . 
Then, a method of carrying out the DFT under such indicated 

15 resolution will be described. 

Setting of the resolution is carried out by the basic 
frequency setting portion 2 represented on the recursive 
discrete Fourier transformation device shown in Fig . 1 and based 
on the resolution information set up here, the Fourier 

20 transformation is carried out such that while the minimum 
frequency of arbitrary frequency band is f 1 and the maximum 
frequency is f2, a desired resolution (frequency interval), 
indicated with any of values which divides an interval between 
fl and f2 to N, is N point. 

25 A method of carrying out the discrete Fourier 

transformation operation with a desiredresolution in a desired 
arbitrary frequency band at every sampling time interval will 
be described in detail. 

In the equations (11 ) and ( 12) , or (19) used for carrying 

30 out the recursive discrete Fourier transformation, cos [2 
7Tk/N] and sin[2 7tk/N] expressed with circular functions 
indicate basic functions relating to analysis point number 
N which should be obtained by operation. 

Therefore, the values of fl, f2 and fs are set up so 

35 that k/N relating to the basic frequency indicated in the 
circular functions coincides with a desired frequency band. 
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Here, it is assumed that basic frequency a of a basic 
function corresponding to the minimum frequency f 1 is a = f 1/f s 
and basic frequency 6 of a basic function corresponding to 
the maximum frequency f2 is 6 = f2/fs. 
5 Then, when the resolution is N points , the basic 

frequency in the frequency band of fl to f2 is as follows. 

(where k = 0 ~ N-l) 



10 Therefore, by replacing the basic frequency k/N in the 
equations (11) and ( 12) or (19) with a following equation (20) , 
it is possible to obtain recursive Fourier transformation 
equation adaptive for transformation with a desired arbitrary 
resolution in an arbitrary frequency band desired to be 

15 analyzed. 

±{(f2-n)k i ... 

fs\ N-l J 

(where 0<£<JV-1) 



Equations for carrying out transformation at an 
20 arbitrary frequency resolution, which corresponds to the 
equations (11), (12) are expressed as equations (21), (22). 



X r {k,t +l)=^X r (k,t) + ^[x{t+N)-x(t)'^ 



fs[ N-l 



25 X,(k,t+l) = X t (k,t) 



"' (21) 
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- (22) 



An equation for carrying out transformation at an 
arbitrary frequency resolution, which corresponds to the 
equation (19), is expressed as an equation (23). 



where A is a positive constant for providing [x(t+N) 
- x(t)] with an amplitude value, and 



Fig. 5 shows the structure for carrying out the Fourier 
transformation at an arbitrary frequency resolution. 

In the same Figure , sampling data sampled at every equal 
time interval is supplied to a Fourier operating portion 6, 
which is a linearly invariable system, in succession, so that 
for example, data x(t+N) is supplied at time t+N. 

At this time, a basic frequency setting portion 7 sets 
up a lower limit frequency f 1 and a upper limit frequency f 2 
of an arbitrary frequency band to be analyzed and then a desired 
resolution N. Then, a basic frequency to be transformed is 
determined according to the aforementioned equation (20) and 
supplied to the Fourier operating portion 6 , which is a linearly 
invariable system. The Fourier operating portion 6 carries 
out IIR filter processing having a transfer function indicated 
by the aforementioned equation (19) according to the basic 
frequency. 

The Fourier operating portion 6 supplies X r (k, t + 1) as 



m 




cos[2jtp]- 7sin[2^?]-z M 
l-2cQs\hrp\~ l +z~ 2 j 



- (23) 




the real part of the Fourier coefficient to N sampling data 
supplied since time t+1 and Xi(k, t+1) as the imaginary part 
thereof. Consequently, the DFT arithmetic operation can be 
carried out in a desired arbitrary frequency band at a desired 
5 resolution at every sampling time interval. 

As described above, recursive discrete Fourier 
transformation operation, in which the resolution is set up 
arbitrarily for analysis, can be executed. This Fourier 
transformation processing can be carried out based on the 

10 equations (11) and (12), or equation (19) for every sampling 
value supplied newly. By changing the value of k to 0, 1, 
2, N-l in that operation, a result of transformation 

processing under each basic frequency set corresponding to 
each k is provided. 

15 Because in the transformation equation for conducting 

recursive discrete Fourier transformation operation, the 
basic frequency set by a predetermined k value is capable of 
obtaining a result of transformation processing independently 
of analysis on the basic frequency set by other k value, a 

20 result of Fourier transformation operation on a specific basic 
frequency relating to arbitrary k value may be obtained 
independently . 

Then, because when attention is paid to a specific 
frequency point or k is assumed to be a specific value, the 

25 value of circular functions in the equations (11) and (12), 
or equation (19) is a constant. Thus, by using the 
aforementioned definition equation (1) for the discrete 
Fourier transformation, which has been employed since before, 
a result of the Fourier operation relating to that basic 

30 function can be obtained more easily than the Fourier 
transformation relating to a specific basic function. 

Because the arithmetic operation for the signal 
component relating to the specific basic function is more 
difficult in case of the FFT which uses butterfly operation, 

35 the recursive discrete Fourier transformation method shown 
here is expected to be widely applied as a signal analysis 
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method, which will be employed in more diversified ways in 
future . 

The recursive discrete Fourier transformation method 
has been described above . Because this transformation method 
5 is capable of selecting a basic frequency preferable for 
analysis, the number of points to be handled is not restricted 
to power of 2 unlike in case of the FFT, so that any positive 
integer can be adopted. 

Although in the above description, it is assumed that 
10 the sampling data to be supplied is real number data, even 
if the sampling data is complex data like shown below, it can 
be expanded . 

Specifically, in case where the sampling data is complex 
data, an equation corresponding to the equations (3) , (4) is 
15 as follows. 

X F (k,t + 1) = |z r (*,/) + l[x r (t + N) + x r (f )]j cos |2^J 

-{zXM^k^^-^Mjsin^] - (24) 

+ |jr r feO + ^k&+^)-x r (0]|sin|2^-J - (25) 

20 (where, A = 4n,0 &k < N -1 ) 

Next, inverse discrete Fourier transformation, which 
is a second embodiment achieved by applying the above discrete 
Fourier transformation technology, will be described below. 
25 Fig. 6 shows an embodiment of the recursive inverse 

discrete Fourier transformation device , to which the recursive 
inverse discrete Fourier transformation method is applied. 



This embodiment will be described with reference to the same 
Figure. It is assumed that the sampling data to be supplied 
is complex data. 

This recursive inverse discrete Fourier transformation 
5 device comprises a data updating portion 11 , which is supplied 
with data sampled at every predetermined time interval, for 
storing supplied new N data (N: positive integer) temporarily, 
a basic frequency setting portion 12 for setting the basic 
frequency for carrying out the inverse discrete Fourier 

10 transformation, a recursive inverse discrete Fourier 
transformation (IDFT) operating portion 13 and amemory portion 
14 for storing computed data temporarily. 

An operation of the recursive inverse discrete Fourier 
transformation device having such a structure will be 

15 described. 

First, supplied complex data is sampled by a sampling 
circuit (not shown) at every predetermined time interval, so 
that discrete data of the real part and imaginary part , sampled 
and quant itized, are supplied to the data updating portion 

20 11. 

The sampling circuit samples data values supplied at 
times t, t+1, t+2, t+3, t+N-1, t+N (N: natural number) 

each having an equal interval so as to generate data stream 
x r (t) + jXi, x r (t+l) + jxi(t+l), x r (t+2) + jXi(t+2), x r (t+3) 

25 + jXi(t + 3) , . . . ,-x r (t+N-l) + jXi(t+N-l) , x r (t+N) + jXi(t+N) with 
data value supplied at each time as sampling value 
corresponding to each time. 

In this manner, the sampling data x r (t) + jxi(t) sampled 
at time t, x r (t+l) + jXi(t+l) at time t+1 -x r (t+N-l) + 

30 jxi( t+N-1) at time t+N-1 -x r ( t+N-1) + jXi( t+N-1) at time 

t+N-1, x r (t+N) + jXi(t+N) at time t+N are supplied to 

the data updating portion 11. The data updating portion 11 
updates complex data of new N data (N: positive integer) of 
the supplied data and stores them therein temporarily. 

35 The N data stored temporarily is supplied to the 

recursive IDFT operating portion 13 and a last FFT operation 
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result stored in the memory 14 temporarily is supplied to the 
recursive IDFT operating portion 13 as recursive data. Then, 
the recursive IDFT operating portion 13 carries out recursive 
discrete inverse Fourier transformation according to 
frequency resolution information set up by the basic frequency 
setting portion 12 in a manner described later and outputs 
a result of its arithmetic operation. 

Then, the operation result subjected to inverse discrete 
Fourier transformation by operation by the recursive IDFT 
operating portion 13 is supplied to the memory portion 14 and 
stored therein temporarily . The temporarily stored operation 
result is supplied to the recursive IDFT operating portion 
13 . Here, the IDFT operation for next supplied data is carried 
out, so that the result of the IDFT operation is used 
recursively. 

Next, the IDFT operation, which is carried out 
recursively, will be described in detail. 

Inverse Fourier transformation X(k, t) at arbitrary time 
t is defined according to a following equation (26). 

X(i ( <)--^"|' X ( n )exp|j2(n-()^J - (26) 

- Wt* ' Wcos [ 2( " < ] -wt** l 2[n - t] f}} 

= X r {k,t)+ jx f (Kt) 

for k = 0,l » N-l 

Data stream x(n) supplied here is complex number, which 
is a number expressed like x(n) = x r (n) + jxi(n) . 

The recursive inverse discrete Fourier transformation 



is carried out for N complex data supplied in this manner and 
its arithmetic operation is carried out through the same 
processing as the aforementioned recursive DFT transformation 
method. 

5 In this operation, complex Fourier coefficient X r (k, 

t + 1) + jXi(k, t+1) at sampling time t + 1 can be obtained according 
to following equations (27), (28) by using a difference between 
newly complex data x r (t+N) + jXi(t+N) sampled and fetched at 
sampling time t+1 and the oldest complex data x r (t) + jxi(t) 
10 of N-data stream and complex Fourier coefficient X r (k, t) + 
jX ± (k, t) already obtained at sampling time t. 




- (27) 



- (28) 



(where, B = -Jn, 0 < N -1) 



The recursive IDFT transformation is carried out in the 
above way. Like the Fourier transformation, the inverse 
20 Fourier transformation is also capable of transforming the 
newest N-point data stream containing new data, which is 
sampled successively at every sampling time interval shown 
in Fig. 2. 

Meanwhile, in case where the sampling data is complex 
25 data, an equation corresponding to the equations (21), (22) 
is as follows. 
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Z r (*,r + l)-j^M + ^('^)-^(')])KC0sj2|[(^^ + /l]j 

-{z,M + ~k(<^)^,(<)]}>< S m{2^[^^^ + /l]}--- (29) 

* 1 (t >t+ i)-{ir I M+^kfr+^)-*.W}*««{2f 

+ {x,(*,0 + ik(< + ^)-^W]}xsm|2i[to^ + / 1 ]j--- (30) 

5 Corresponding to the equations (29). (30), an equation 

for carrying out inverse transformation at an arbitrary 
frequency resolution is expressed in the form of equations 
(31), (32). 

AT r (*,f + l)-{*>,0H^ 
10 ^kO + lfcfr^ (31) 

-{z r M4M' + ")-^ (32) 

It should be understood that many modifications and 
15 adaptations of the invention will become apparent to those 
skilled in the art and it is intended to encompass such obvious 
modifications and changes in the scope of the claims appended 
hereto . 
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